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We present a new chiral quark model for mesons involving a heavy and a light 
(anti-) quark. The model relates various combinations of a quark - meson coupling 

predicted and some used as input. The extension from other similar models is that 
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the present model includes the lowest order gluon condensate of the order (300 



O 

I MeV) determined by the mass splitting of the and the 1 heavy meson states. 



£NJ ■ Within the model, we find a reasonable description of parameters such as the decay 

> 

constants Jb and fjj, the Isgur-Wise function and the axial vector coupling gA in 



chiral perturbation theory for light and heavy mesons. 



Qh; I. INTRODUCTION 

Oh! 

While the short distance (SD) effects in hadronic physics are well understood within 



perturbative quantum chromodynamics (pQCD), long distance (LD) effects have been hard 
to pin down. Lattice gauge theory should be able to solve the problem, but the calculations 
are often very difficult to perform. QCD sum rules might also give the answer in some cases. 
Still, in some cases it has been fruitful to use various QCD inspired models and assumptions. 
In the light quark sector, low energy quantities have been studied in terms of the (extended) 
Nambu-Jona-Lasinio model (NJL)|l[] and the chiral quark model (xQM)||, which is the 
mean field approximation of NJL. 

In the xQM, the light quarks (u, d, s) couple to the would be Goldstone octet mesons 
(K, 7r, rj) in a chiral invariant way, such that all effects are in principle calculable in terms of 
physical quantities and a few model dependent parameters, namely the quark condensate, 
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the gluon condensate and the constituent quark mass 0, |, More specific, one may 
calculate the coupling constants of chiral Lagrangians in this way. 

In this paper we will extend the ideas of the chiral quark model of the pure light sector 
]2|, |3|, |], [|] to the sector involving a heavy quark (c or b) and thereby a heavy meson. Such 
ideas have already been presented in previous papers |6|, |7|, §|, 0- Also in this case, one may 
calculate the parameters of chiral Lagrangian terms, where the description of heavy mesons 
are in accordance with heavy quark effective field theory (HQEFT) . In the present paper we 
will extend the ideas of H 0, H, [§ to include gluon (vacuum) condensates. The motivation 
for the inclusion of gluon condensates is that this works well in order to understand the 
AI = 1/2 rule for K — > 2n within the xQM @, |5| and within generalized factorization [10 



Furthermore, it allows us to consider effects related to the gluonic aspect of rj' as consid- 
ered in some aspects of D- meson decays |p!2| , and also to calculate gluon condensate 
contributions to B — fi-mixing ]13| . 



Having established our heavy - light chiral quark model (HL^QM), we can integrate out 
the light and heavy quarks and obtain chiral Lagrangians involving light and heavy mesons 



n, 15]. Chiral perturbation theory (%PT) based on such Lagrangians works in the pure 



strong sector. In order to define the model and its parameters, we have to integrate out the 
quarks, and we will find some typical divergent loop integrals. We will relate all divergent 
loop integrals to some physical parameters, as was done in [El. This means in particular 
that we will treat quadratic, linear, and logarithmic divergences as different. If we need to 
calculate a divergent integral we will do so in dimensional regularization, although various 
regularization procedures might be used. Thus, the regularizing prescription for divergent 
diagrams is to be regarded as a part of the model. Still, even if integrals are divergent, the 
effective UV cut-off scale is, as for xPT, the chiral symmetry breaking scale A x ~ 1 GeV, 
where also the matching of pQCD and HL^QM is performed. This is also considered a part 
of our model construction. 

The paper is organized as follows: In the next section (O) we describe the HLxQM, 



and in section |T| we consider bosonization in the strong sector and of the weak current 
respectively. In section [IV] we discuss the relations between physical and model dependent 
parameters. In section [VJ we discuss the Isgur-Wise function, and in section [VT] we bosonize 
the 1/rriQ terms. Section VII contains the presentation of some necessary chiral corrections 
for our numerical analysis. Finally in section |V1I| we discuss our results. Loop integrals are 
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listed in Appendix [A]. In Appendix [B] we list some transformation properties of the involved 
fields. 

II. THE HEAVY - LIGHT CHIRAL QUARK MODEL (HLxQM) 

Our starting point is the following Lagrangian containing both quark and meson fields: 



C = C 



HQEFT 



Int 



(1) 



where 16 



£hqeft = Q v iv ■ DQ V + 



2m Q V V 2 



-2\ 



D ± Y eS ) Q v + 0{m- Q 



(2) 



is the Lagrangian for heavy quark effective field theory (HQEFT). The heavy quark field 
Q v annihilates a heavy quark with velocity v and mass uiq. Moreover, is the covariant 
derivative containing the gluon field (eventually also the photon field), and a-G = 'G a ^ v t a , 
where = i[ r y IJ ', r y v ]/2, G a is the gluonic field tensors, and t a are the colour matrices. 
This chro mo-magnetic term has a factor Cm, being one at tree level, but slightly modified 
by perturbative QCD.(When the covariant derivative also contains the photon field, there 
is also a corresponding magnetic term ~ a ■ F, where F^ v is the electromagnetic tensor). 
Furthermore, (iD^)^ = C D (iD) 2 — C K iiv ■ D) 2 . At tree level, C D = C K = 1. Here, C D is 
not modified by perturbative QCD, while Ck is different from one due to perturbative QCD 
corrections [[T7J . 

The light quark sector is described by the chiral quark model (xQM), having a standard 
QCD term and a term describing interactions between quarks and (Goldstone) mesons: 



xQM 



qii^D^ - M q )q - m{q R Y)q L + q L Y,q R ) 



(3) 



where q T = (u, d, s) are the light quark fields. The left- and right-handed projections qi and 
qR are transforming after SU(3)l and SU(3)r respectively. M. q is the the current quark 
mass matrix, m is the (SU(3) - invariant) constituent quark mass for light quarks, and 
X = exp(i ^2 j is a 3 by 3 matrix containing the (would be) Goldstone octet (it, K, rj) 



£ = e in// where n = — a (x) = 

2 v 2 



7t" I m 
K~ 



71 



21— _l_ M. tf° 
~V2 ~r V6 



T^ 8 



(4) 
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The xQM has a "rotated version" with flavour rotated quark fields x given by: 

Xl = ; xr = €qr ; f • f = s • ( 5 ) 

In the rotated version, the chiral interactions are rotated into the kinetic term while the 
interaction term proportional to m in @ become a pure (constituent) mass term || || : 

£xQM = X [Y(iD^ + V M + 75 AJ "Hx- , (6) 

where the vector and axial vector fields V M and are given by: 

v^e^e+e^) ; ^= - - e^e f ) , (7) 

and M q defines the rotated version of the current mass term: 

M q =M v q + Mf lh , where (8) 

M? = l@Mtf+ tM£) and Mf=-l$M$-£M£). (9) 

Here L is the left - handed projector in Dirac space, L = (1— 7s)/2 and R is the corresponding 
right - handed projector, R — (1 + 7 5 )/2 . The Lagrangian (§) is manifest invariant under 
the unbroken symmetry SU(3)y. In the light - sector, the various pieces of the strong 
Lagrangian can be obtained by integrating out the constituent quark fields \i an d these 
pieces can be written in terms of the fields A^ , and which are manifest invariant 
under local SU(3)v transformations. For instance, the standard 0(p 2 ) kinetic term may 
(up to a constant) be written as Tr [A^ A u ]. This is easily seen by using the relations 

A = ~y (D,tf) e = + (d^) £t , (io) 

where is a covariant derivative containing the photon field. The vector field V M transforms 
as a gauge field under local SU(3), and can only appear in combination with a derivative as 
a covariant derivative (id 1 * + V M ). 

In the heavy - light case, the generalization of the meson - quark interactions in the pure 
light sector xQM is given by the following 577(3) invariant Lagrangian: 

&rU = -G H [XaHZQv +Q~vH:Xa] + \W* H«] , (11) 

where Gh and G3 are coupling constants, and H% is the heavy meson field containing a spin 
zero and spin one boson: 

H a v =p+(p; 7 m - iPti*) , 

h° = 7°(^)V = [(p;) V - kpsH] p + , (12) 
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FIG. 1: Feynmanrule for the light quark -soft gluon vertex. 



where 

P± = {\±i>)/2. (13) 

are projection operators. The index a runs over the light quark flavours u, d, s, and the 
projection operators have the property 

P±1^P± = ± P± P± ■ (14) 

Note that in || 0, || ||, Gh = 1 is used. However, in that case one used a renormalization 
factor for the heavy meson fields H v , which is equivalent. 

The fields P 5 (P At ) annihilates a heavy-light meson, CT(1 - ), with velocity v. The interac- 
tion term in the Lagrangian (equation (|ll])) can, as for the xQM, be obtained from an NJL 
model. In the NJL model one starts with the Lagrangian for free quarks and four quark 
operators, thought to be generated by gluon exchange(s). Taking the heavy quark limit for 
heavy quarks one can obtain (after some manipulation) the interaction term ( |TT1) from the 
four quark operator . This has been done in (-as for the light sector Jl|). 

In our model, the hard gluons are thought to be integrated out and we are left with 
soft gluonic degrees of freedom. These gluons can be described using the external field 
technique, and their effect will be parameterized by vacuum expectation values, i.e. the 
gluon condensate (^-G 2 ). Gluon condensates with higher dimension could also be included, 
but we truncate the expansion by keeping only the condensate with lowest dimension. 

When calculating the soft gluon effects in terms of the gluon condensate, we follow the 
prescription given in The calculation is easily carried out in the Fock - Schwinger gauge. 
In this gauge one can expand the gluon field as : 

A^k) = -^G p ,(0)^(k) + .... (15) 

Since each vertex in a Feynman diagram is accomplished with an integration we get the 
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Feynman rule given in figure |I[ The gluon condensate is obtained by the replacement 



We observe that soft gluons coupling to a heavy quark is suppressed by 1/tuq, since to 
leading order the vertex is proportional to v fl v u G af11 ' = 0, t> M being the heavy quark velocity. 

III. BOSONIZATION WITHIN THE HLxQM 



The interaction term Cj nt in (|TT|) can now be used to bosonize the model, i.e. integrate 
out the quark fields. This can be done in the path integral formalism and the result is 
formally a functional determinant. This determinant can be expanded in terms of Feynman 
diagrams, by attaching the external fields H%, H%,V^, and M^' A of section II to quark 
loops. Some of the loop integrals will be divergent and have to, as for the pure light sector 
|2|, Q ||, to be related to physical parameters. The strong chiral Lagrangian has the 
following form fTJ, |T5|, fZOj, gTJ, 0, jEJ : 



Cstr = -(1 + — )Tr \H~ a (iv ■ D)H a ] + (A Q + -^-)Tr \JT a H a ] 

m Q - m Q 

+ (1 + ^-)Tr [H a H b vXa] - (9A ~ —)Tr [H a H bl , l5 Al] 



+2X 1 Tr 



H a H b {M^) ba + 2X[Tr [H a H a ] (M^) bb + 



A2 -Tr \7T a cr^H a a al3 ] + ^-Tr [H a o af3 iv ■ DH a a aP ] 



ArriQ m Q 

- ^-Tr [~H a <T a %Vg a a af3 H b ] + ^-Tr \W allll ^ a H h } + .... (17) 



where the ellipses indicate other terms (of higher order, say), and contains the photon 
field. The 1/mq terms will be discarded in this section, but will be considered later in 
section VI. Note that the term proportional to Aq = Mh — tuq is absent in most articles 
considering the mesonic aspects only. The quantities Aq, Ai, A2 and 5q are related to the 
masses of the heavy mesons. The trace runs over the gamma matrices (only). Note that in 
some conventions there is an extra factor Mh at the right- hand side in (|17]) (our notation 
is the same as used in || [7|, || [14], |15[ ) . 

The Feynman diagrams responsible for the kinetic term and the mass difference (self 
energy) term ~ Aq in the meson Lagrangian is shown in figure 0. A calculation of these 
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two diagrams (at zero external heavy meson momentum), leads to the identification 

-iG 2 H N c (-h + mI 3/2 + ^ G 2 )) -±- = -A Q , (18) 

where we have added the last term of (|TT|). Keeping the heavy meson momentum to first 
order (subsequently to be interpreted as the derivative of the heavy meson field) , we obtain 
the identification for the kinetic term: 

-iG 2 H N c (7 3/2 + 2ml 2 + ^r(—G 2 )) = 1 . (19) 

l\ c IT 

We have denoted the divergent integrals by / and the finite integrals by k. They are 
defined in Appendix [A|. The quadratic - , linear - and logarithmic - divergent integrals are 
denoted Ji, I 3 / 2 , and l 2 respectively. 

The relation flT9| ) is also obtained by comparison of the loop integral for diagram in figure 
|3] with the vector field attached to the light quark. This must be so because of the 
relevant Ward identity, but it is also realized by explicit calculation. (Note that the total 
covariant derivative is iD^ = iD^ + V M in the quark sector and iD^ = iD^ — in the meson 
sector since the meson fields transforms as an anti triplet under SU(3)y. See Appendix |B] 
for details) 

From the same diagram, with the axial field An attached, we obtain the following iden- 
tification for the axial vector coupling : 

g A =iG 2 H N c (\l 3/2 - 2ml 2 + i«o - -^-(— G 2 )) , (20) 

6 6 JM C IT 

and similarly when attaching MY 

2X 1 =iG 2 H N c (I 3/2 - 2ml 2 + 2k - ^(-G 2 )) , (21) 

JM C 7T 

Within the full theory (SM) at quark level, the weak current is : 

Jf = qJ Ll a Q (22) 
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V a ,A a 



FIG. 3: Coupling to vector and axial vector current 



where Q is the heavy quark field in the full theory. Within HQEFT this current will, below 
the renormalization scale \i = irtq (= mj, m c ), be modified in the following way: 



Jf = Xa£ f T a Qv + 0{m Q l ) 



(23) 



where |T6] 



T a = C 7 (/i) 7° L + CM v a R . 



(24) 



The coefficients C 7) „(/^) are determined by QCD renormalization for fi < tuq. They have 



been calculated to NLO and the result is the same in MS and MS scheme p4|. Corrections 
to the weak current of order 1/mq will be discussed in section |VT] . 

The operator in equation (^) can be bosonized by calculating the Feynman diagrams 
shown in figure |] : 

j;(Bos) = Jf(0) + Jf (!) + ••• (25) 
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FIG. 4: Diagrams for bosonization of the left handed quark current 



The currents Jf(0) and J" (1) correspond to zero and one axial field attached to the loop, 
the dots represents terms with more boson fields and gluon condensates. We obtain to zero 
order in the axial field and first order in the gluon condensate: 



Off, 



J°(0) = ^Tr 



where 



k 2 ,ol 



a H = - 2iG H N c -h + mh/2 + ^(—G 2 ) 



iV c N vr 



(26) 



(27) 



Note that this expression is also related to I/G3 in eq. (jig) . However, in the present case 
only Gh in first power is involved. 

To next order in the chiral expansion we obtain the current 



JJ(1) 



where the quantities a H ^ are given by 



-Tr 
2 



(28) 



a: 



(1). 

#7- 



ft 



(1). 



2iC H iV c (~I 3/2 + ^0 - 2mJ 2 + ^-(^ G 2 >) 



(29) 
(30) 



We observe that these quantities, as fu in (|33|) , contains only one power of the coupling 
Gh, in contrast to the strong sector. 

The physical meaning of the ft#'s becomes more clear by considering their contributions 
to the physical quantities /#*, and the semileptonic form factors f±(q 2 ). The coupling 
an in ( p6|) is related to the physical decay constants /# and /#*, in the following way (for 
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= =X 



FIG. 5: Feynman diagrams for the H — > ir transitions, both J" (0) and J° (1) defined in (^) and 
(|28|) respectively, contributes to the diagram to the left 



H = B,D): 



(0|«7 Q 75&|#> = -2 <0|J^|iT> = iMnfnv , (31) 
(0|M 7 a &|#>= 2(0|J a Q |tf> =M H ,f H ,e a . (32) 



Taking the trace over the gamma matrices in (EBH, we obtain a relation for an and the 
relations between the heavy meson decay constants fjj and /#* (for H = B,D) : 

aH c 7 (^)+a(/i) c» ' 1 j 

where the model dictates us to put // = A x . The form factors f+(q 2 ) and f-(q 2 ) are defined 
as : 

(7r + G^ 7 a (l-75)^> = 2(7r + ( P7T )\Jf\H) = U(q 2 )(p H +Pn) a + f-(Q 2 )(PH-PnT (34) 

where p H = Mnv a and the index a corresponds to quark flavour u and g M = p H — k%. 

The diagrams in figure |5| contributes to the H — > 7r transition. J?(0) and is 
responsible for the diagram to the left and the diagram to the right involve the strong 
Lagrangian term ~ in ([T7|). A calculation of the diagrams gives : 

U(q 2 ) + f-tf) = 7=^ (C, + C V - g A C 7 ) a H , (35) 
2n f I 2\ . n ^/Mh ( 9A&H (1) 



/+M -/-M = -^^ + ^J, (36) 

where we have neglected terms of first order in v ■ (where a H \ contributes). This means 
that the equations are only valid near the "no-recoil point" , where v -k n — > and q 2 — > g^ ax = 
(M 2 H + ml). From equation © and @ we see that (f+(q 2 ) + f-(q 2 ))/(f+(q 2 ) - f-(q 2 )) ~ 
1/Mh, which is the well known Isgur-Wise scaling law [^jj. (The 1/v ■ k n term in (^) is 
due to the if* pole). 
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IV. CONSTRAINING THE PARAMETERS OF THE HLxQM 

Within the pure light sector the quadratic and logarithmic divergent integrals are related 
to f n and the quark condensate in the following way || [|, [|, p6| : 

f* = -^NJ 2 + 1 ^(^), (37) 



(qq) = -AimN c I x - -J-(^G 2 > , (38) 

12m 7r 



which are obtained by relating loop diagrams to physical quantities as for the eqs. ( |T9| ) and 
(p0|). (Here the a priori divergent integrals have to be interpreted as the regularized 
ones) As the pure light sector is a part of our model, we have to keep these relations in the 
heavy light case studied here. In addition, in the heavy-light sector the linearly divergent 
integral I3/2 will also appear. As 1\ and I 2 are related to the quark condensate and f n 
respectively, the (formally) linearly divergent integral J3/2 is related to 5g^ = 1 — g^, which 
is found by eliminating I 2 from eqs. (19) and ( p0|) : 



$9A = ~\iG 2 H N c (l 3/2 + no) ■ (39) 

Note that the gluon condensate drops out here. Within a primitive cut-off regularization, 
J3/2 is (in the leading approximation) proportional to the cut-off in first power ||. Within 
dimensional regularization it is finite. We will keep 1 3/2 as a free parameter to be determined 
by the physical value of gj±. 

Eliminating ^3/2 from the eqs. ( |H)D and fl2"D| ) and inserting the expression for 1 2 obtained 
from (|37|) we find the following expression for Gh '■ 

G\ = m [ l N + *g A ) where ^_ 

9 f 2 1 rn 2 N c _ JTWOs /-<2\ 32 

Z ^7T 47T m 2 \ 7T ^ / 

Note that Gh has dimension (mass) -1 / 2 . 

In order to constrain the parameters further, we will consider the parameters Ai, A2 and 
5q related to the meson masses. The gluon condensate can be related to the chromomagnetic 
interaction : 

^ g{h) = 2M^ CM{H ^4 a ' GQvlH) ' (41) 

where the coefficient Cm(^) contains the short distance effects down to the scale /i and has 



been calculated to next to leading order (NLO) pTL p8| , and can be found in table H (2M H 
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is a normalization factor). The chromomagnetic operator is responsible for the splitting 
between the 1~ and 0~ state, and is known from spectroscopy, 

H%(H) = 3A 2 = ^m Q (M H * - M H ) . (42) 

An explicit calculation of the matrix element in equation (|4*T|) gives 

fJ'G = ?72 — (— G 2 ) , where m = ^^C M (A x ) . (43) 

Combining eq. ( ffOD and eq. ( f4"3|) we get the following relations : 



^ G) "^Tp' G ""7f p ' (44) 



where the quantity p is of order one and given by 

;i + 3<u) 



v , 2m2 

4(l + ff) 

In the limit where only the leading logarithmic integral I2 is kept in (|19"D, we obtain: 

t(0 ) \/2m 



<M - 1 , p - 1 , G H - G™ = ^— . (46) 

Jit 



Note that £4 = 1 is the non-relativistic value [15 



The quantity 5q is found from the kinetic term as: 

6q = ^l^( H \Q~v(D ± ) 2 cS Q v \H) , where (D ± f ee =D 2 - C K (v ■ D) 2 (47) 

This quantity can easily be calculated in our model, and a direct calculation gives the 
expression: 

/4 = l G 2 H N c m 2 {h /2 + 2mI 2 + ^G 2 ))+r 1 ^{^G 2 ) 

+-C K G 2 H { -MmN c h - — (-G 2 ) 1 , where m = — + — . (48) 
4 \ 12 N 7T 7 J /d 48 64 K ' 

Note that the (v ■ d) 2 part of the kinetic term cancels the two heavy quark propagators, 
and the light quark condensate (|3"%| ) appears naturally Eliminating the divergent integrals, 
using equation (0), ( |3"8"D and (|4"3"D, some of the G^'s can be included in physical parameters, 
and we obtain 

5q = /4 = -/4 -m 2 --C K {qq)G 2 H . (49) 
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Note that the last term (originating from (v ■ D) 2 Q V at quark level) gives a vanishing 
contribution in the free quark limit when Gh —* 0. Using the expressions for (^-G 2 ) and 
G 2 H in (§4]), we find an expression for g A in terms of m, fj, 2 and (qq). However, as fj, 2 is not 
very well known, and Ck is known only to leading logarithmic approximation, we will not 
try to determine $a by this relation. 

From eqs (|2T]), (|37|) and (g|), we find 



2Ai = -(3&i-i; 



(97T - 16)//| 



(50) 



384?72m 2 

In the limit (fi6| ) we obtain 2Ai — >• 1, as expected. The parameter Ai is related to the mass 
difference M Hs — M Hd . Unfortunately, this cannot be used at the present stage to constrain 



the parameters in our model because this quantity has large chiral corrections |23 
Using equation fllTf ) and (^) we can write an as 

G H ( (qq) ~ 1* U - 2 



an 



p ~^ 16m 2 



or equivalently, using the eqs. (|38|) , (|39|) and (^Tj) as 



GW = - m 1 - g A ) + \ ' — + G 
2 192 1] 2 m 



H 



(-G 2 ) 

IT 



( qq ) m 2 
2m + 8vr 



(51) 



(52) 



Note that the relation ( p3[ ) gives fn > fn*, which is not correct experimentally. Adding 



l/rriQ corrections we correctly reproduce fn* > fn pi 



Combining ( |33"D with ([51]), we obtain |12j in the leading limit (taking into account the 
logarithmic and quadratic divergent integrals only, and let C 7 — ► 1, C v — *■ and £4 — > 1 as 



m 



UV^m 

which gives the scale for /# (It is, however, numerically a factor 2 off for the 5- meson). 



(53) 



Using the relations in equation (0), ( p0|) and fl3"T|) we obtain for a^, and 



(1) 



a 



(i) 

(1) 



2^4 

Gh 

A G ( l 1 



mAf, 



256m 3 7r 



(54) 
(55) 



where the latter may also be written: 



yAJ 3 H 4vr 



mN r (tt 
+ 



256?72m 2 



(56) 
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FIG. 6: Loop diagrams for bosonizing the b — * c current, = 7^, = 7^75 



Combining equation fl35|) and (|36|), we find (up to terms of first order in v ■ k n ): 



9AfH-y/M H M H . „ g A VM^ 1 r [Mh~* 



Utf) = - *~V, , -T~nnr{fH-9AfH*\H7-} . ( 57 ) 



where we have used equation fl33|) and 



V. THE ISGUR-WISE FUNCTION 



The Isgur-Wise function relates all the form factors describing the processes 

B — > D(D*) in the heavy quark limit : 



yj Mb Mb =^)(v» + v ),. = vv 



(D*\Q cv riMB)_ = ^ ^ ^ _ e , M(1 + w)] ; (5g) 



= a^ uXa e* u v' x v a , (59) 



where Q c „ and Qb v are the c and b quark fields within HQEFT. The Isgur-Wise function 
(IW) can be calculated straight forward by calculating the diagrams shown in figure (^j). 
The result is 

= -A- +,,.(«), (60) 
l + u 

where p is given in fl4"5p and r(u;) is the same function appearing in loop calculations of the 
anomalous dimension in HQEFT : 

r(u) = A_ In (w + vV - l) (61) 

We see that £ is normalized to 1 at zero recoil. Note that equation (|60| ) is only valid in the 
limit of equal gluon condensate and the coupling Gh in the B and D sector. The derivative 
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of the IW function at zero recoil is : 



dw 



w=l 



2 + r 



(62) 



In the limit (pE6|), the IW function takes the simple form 

£(u) = r(u) and 



dw 



w=l 



1 

3 



(63) 



Numerically, the full result ( |60|) is only a few percent away from (|63|) . Adding the short 
distance QCD effects to ( |60|) will slightly modify our result ||27|| . For the IW function 
describing B —>■ D* transition, the derivative at ou — 1 gets a contribution —0.07 from QCD 
corrections such that £'(0) ~ —0.40. 



VI. 1/rriQ CORRECTIONS WITHIN THE HLxQM 

A. Bosonization of the strong sector 

The new terms of order 1/mg in flTTDare a consequence of the chromo magnetic interac- 
tion (the second term in equation (^)), which break heavy quark spin symmetry, and the 
kinetic term (the third term in @). Calculating the diagrams of figure gives the following 
identifications : 

£i = -iG 2 H N c (h - mI 3 /2 + 2m 2 I 2 - uik (64) 

- 2 / 2 ) + ^(-G 2 )) 
4 N c n J 

( 1 4 2 _ 

9i = -iG 2 H N c (--/! + mI 3/2 + -m 2 / 2 + mK - -k (65) 



4 A c 7T 

22 = C M (AJ_^G H (_ G ) ; £2 = -- ■ (66) 

As the 1/rriQ terms break heavy quark spin symmetry, the chiral Lagrangian in (|i~7|) will 
split in H(0~) and H*(l~) terms respectively. However, by allowing for a flavour and spin 
dependent renormalization constants we can write the Lagrangian in the compact form: 

C = -Tr \H*(iv ■ V ba - AqYHI] - g A Tr [HlHfr^Al] , (67) 
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FIG. 7: Diagrams responsible for 1/mQ terms in the chiral Lagrangian 



where zT 5 ^ = i5}, a D^ — V^ a and we have redefined the H fields as H = if r y 'Zh, where Zh 
and the new coupling g A are now defined : 

r-i _ i , ei-2d M e 2 Ar _ v , a4 



where 



Z^ 1 = 1 + -i ^ Aq = Z# I Aq + -^-2- I (68) 

9 A = 9a (l (ei - 2d A e 2 ) I {g x - d A g 2 ) (69) 

V m Q J m Q 



3 for 1 for H*H coupling 

d M ={ d A =l (70) 

•1 for T —1 for H*H* coupling 



The last term in equation (|T7|) gives a splitting in the mass of the H*(l ) and H(0 ) state. 
This term can be absorbed in a redefinition of Aq, namely Aq — > Aq — c!mA#/4, where 
A H =M H * - M H (Note that p? G = 3m Q A H /2). 

Eliminating the divergent integrals in (|65D and flBTf ) using equations (|3~?| ) (|38|) and (p9|), 
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we can rewrite £\ , g\ and g 2 as : 



2 / (??> , ,2 , V , C K (qq) 2 



4m lo7r lb m 



+^|(-G 2 )) where r^-L^ + 8 - 3tt) (71) 
m z 7r / 128 



,. r ,, .'{QQ) | £ | iVcm 2 (37r + 4) C K (qq) 2 , 
* " m - G ^( 12m" + T + 48^ lb" ( ^T + 3/ - 



+ ^ ( «s G2) \ where 1 ((7 _ 2n) (?2) 

m/ 7T / b4 



_ (7T + 4) 

^ 2 " (^T2) bm" • (73) 

The masses of the particles now before SU(3) breaking terms are taken into account 

M(O-) = mQ + A Q + ittlM (74) 

m Q 

M(l-) = m q + A Q+ (3 ^ Q ^ ) . (75) 

B. Bosonization of the weak current 

In HQEFT the weak vector current at order 1/mq is [P|: 

^ = E ^) J " + 2^" E *i W + 2^" T, A MT k a , (7b) 

i=l,2 Q j ® k 

where the first terms are given in (^) and (p4|), the B/s and A^'s are Wilson coefficients, 
and the O^'s are two quark operators 

Of = qil a t^Qv , O? = g L7 Q (-^ • D)Q V , 
2 a = q L v a iIJQ v , 5 a = g^ a (-™ • D)Q V , 

3 « = q L iD*Q v , 0% = q L {-iD a )Q v , (77) 

The operators are nonlocal and is a combination of the leading order currents Jj and a 
term of order 1/mg from the effective Lagrangian @: 

/ ; ' - i [ dy A T{J?(0),O kin (y)} , 



2m c 

^- = i I AT{J 2 Q (0),O kin (l/)} 
^- = i j AT{Jf(0),O mag (y)} 



Q 

J 4 _ • / j4 



2m. 



Q 



1 AT{J 2 Q (0),O mag (y)}, (78) 
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where 



Okin = 77— Qv (iD±) 2 cS Q v , and O mag = - y^—Q v cr ■ GQ V . (79) 
ZTfin tiffin 



Bosonizing the O® operators in (|77 ) we obtain: 



^TriT^H^iah. + alv,)} (80) 



where r Q/i is defined : 



= B x 7 a L7 M + B 2 v a R^ + B 3 g ail R 

+ B 4 -f a Lv^ + B 5 v a Rv^ 1 + B 6 g a ^ R (81) 
The coefficients and ag can be written after the use of ( |3~7| ) ( |38"D and fl55|) : 

a * = 2G "T \ 2^ ( p - 1} + 6i^ ( 7 G > j = J aH + T { qq) (82) 
«S = «2 + = ^a, + |<Mot> (83) 

Bosonization of these non-local operators is straightforward within our model and the 
result is 

(0\(AiT" + A 2 T?)\H) = --0-Tr{^T a H v } (84) 

(o|(a 3 t 3 q + a a t:)\h) = ^MK Tr {^r a H v } , (85) 

where r Q is given in Q2~3]). 

When 1/rriQ terms are included, we find that fl33l) is modified in the following way: 

fu^mi = a H (C, + C v )(l-%^) 

+ BtoZ + b v <4 _ (c +Cv) (AjiM (86) 

2m Q G^rrtQ 



fu* \/~Mu* = a H C y (l- 



ei + 2e 2 , 



+ *x ^ (3^ + ^) 



^7 5^ (87) 



where 



2m Q 3G H m Q 

5 7 = 2Bi — AB 2 — B 3 + B± — B5 — Bq , 
B v = B\ + B 2 + -B3 + B4 + B 5 + Bq , 

5! = - 2Bi + B 3 -B A + B 6 and 5* = B 1 + B A . (88) 
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The ratio between the coupling constants is : 



ju C 7 + C v [ rriQ \ 3a H G Hi 

+ 2n^ ^ " B>1 + {B: ~ Bv)a ^ ■ (89) 

VII. CHIRAL CORRECTIONS 

Chiral corrections are numerically comparable with the 1/toq corrections. While the 
relevant mass scale of I/tuq corrections are Ah = Mh* — Mh (cfr. eq (^)), the relevant 
mass scale of chiral corrections are 5 = Mh s — Mh u d - Clearly Ah ~ 5 numerically. In this 
section we will consider just the chiral corrections to fn and g*[ Hn which are necessary to 
include in the numerical analysis. These corrections have been calculated by many groups 



21, 30, 31 1 , the result is in MS scheme, including counter terms : 



fk^M HuM = f H: ^M H:d = 

«^(l-^^^{-<(l+^)+<(ln<//i 2 + ^ln^)(l + 3^)|>) , (90) 



flVM^=f H: JM H: 



®h[1- [-^(1 + 9a) + m 2 K (\nm 2 K /fj, 2 + A In i)(l + 3<&)j 



32vr 2 / 2 9 



-nil 



(91) 



Note that this result is independent of Aq. We have ignored terms proportional to m 2 and 
used the mass relation m 2 g = Amj^/S. The counter terms needed to make the expression 
above finite originates from the following terms in the weak current : 

MM) = yTr[&r^ c Aia + ^Tr[^H vb )Mt c , (92) 

where T M is given in (0). Moreover, uoi = —AXi/Gh, where Ai is given in equation (pT|). 
We ignore u[ since it is subleading in 1/N C . (This is similar to what is found for 2Li — L2, 
L4 and Lq in the light sector ||). Therefore, in the limit where we neglect u, d quark masses 
f H does not depend on counterterms. This also happens for g~[ Hn . The chiral corrections 
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to g_x can be calculated with the formula listed in appendix [X| : 

X f 1 2 1 ( 35 2 , m K , 8 2 , 4 5 > 



An^f -2FK 8 /A Q )^))} . (93) 

-F(x) is defined in ( |A15| ). This result coincides with the one in for Aq = 0. It turns out 
that for Aq ~ (0 — 0.5) GeV, the chiral corrections vary with less than 1%. Therefore we 
will simply ignore Aq. 

Note that in the expressions for the chiral corrections to fn and Qa considered in this 
section, the I/tuq corrections of the preceding section are not included. However, both chiral 
and 1/uiq corrections are of course included in our numerical analysis in the next section. 



VIII. NUMERICAL RESULTS 



As we have seen in the previous sections, our bosonizing procedure puts restrictions, in 
the form of relations between the the model dependent parameters m, Gjj, (^-G 2 ), (qq) 
and the measurable parameters (quantities) such as f n , g^, [1q, and fi 2 . We will use a 
standard value of the quark condensate, (qq) = —(0.240 GeV) 3 . (It is not clear if the quark 
and the gluon condensates defined in this paper are the same as those in QCD sume rules). 
In principle, we have enough relations to fix the model dependent parameters. However, 
some physical quantities are still relatively uncertain, which means that the values of the 
model dependent parameters cannot be given a very precise value. 

In principle, Aq = Mfj — mQ is also a parameter of the model, which enters if we calculate 
diagrams with external momenta ( v^Aq will then be a part of an external momentum). 
However, the way we are bosonizing here, the external fields (V, A and H v ) carry zero 
external momenta. Then Aq will not enter our loop integrals, and eq. (113) is irrelevant (so 
far) within our model. In the case of chiral corrections Aq could play a role. However, as 
we have seen in the previous section, Aq does not enter in the chiral corrections to fjj and 
plays a very little numerical role for the corrections to g^. 

In the chiral corrections to fn and g^, we have consequently used MS as in ||. In pure 
xPT the U (MS + 1)" scheme is used. We have explicitly checked that the numbers in table 
fl| and |TTT| can be reproduced in U (MS +1)" with a small change in m and g^ (the bare 
coupling constant). 
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The weak decay constants of heavy mesons have been calculated by many groups. Typical 
results from lattice calculations are f B = (200 ± 30) MeV and f D = (225 ± 30) MeV 
H. QCD sum rules gives f B = (180 ± 30) MeV and f D = (190 ± 20) MeV @ and 



NRQCD gives f B = (147(ll)i 12 (9)(6) MeV||34j|. In order to constrain our parameters we 
will use the following combinations from QCD sum rules which have been evaluated rather 
accurately! 



fBfB-y/2gT Bn M B /f„ = (0.64 ± 0.06) GeV 2 , 

f D fD*V2g°* D *M D /U = (0.51 ± 0.05)GeV 2 , (94) 

where g^* HlT (for H = B,D) is the chiral coupling with chiral and l/rrig corrections 
included. The left-hand side of (0) - for H = B,D respectively- is a function of m and 
(the uncorrected) g^. Therefore fl9"I|) gives as a function of the constituent light quark 
mass m, which has been plotted in figure [3] and [| Thus, we can then plot f B , fo and other 
quantities as a function of the light quark mass. We may further use explicit values for f B 
and fo to determine a value for m in the B- and D- meson sector separately. However, 
because we consider the Isgur-Wise function which involve both B- and D- mesons, we need 
a unique value of m, which give a reasonable value of f B and fo simultaneously. As can be 



seen from the plot in figure |10|, this can be accomplished by taking 

m = (220 ± 30) MeV . (95) 

This is consistent with the value m = (200 ± 5) MeV used in the pure light sector in 
order to fit the AI = 1/2 rule for K — > 2ir decays. 

In table |I| and table [TTT] we have listed some of the predictions of HLxQM for decay 
constants and counterterms. The input parameters have been listed in table |. It could be 
argued that the bare parameters, g^, Gh, {^G 2 ), listed in table [TJ and table fTJ should be 
equal. As statet earlier in this model m has to have an unique walue in both the D- and 



.B-sector. In order to fit our model to the result from QCD sumrules (equation (94[)), we 



have to alow for a different walue of g^, Gh, (^-G 2 ) in the B- and -D-sector. From table |II 



and table [TT| we see that these parameters agrees within errorbars in the two sectors. 

The inclusion of the counterterm for f Bs is crucial, putting to\ = would give /b s //b = 
1.29 ± 0.03 and /d s //d — 1-34 ± 0.04 which is much higher than most lattice estimates |36| 
f B Jf B ~ f Ds /f D ~ 1.15 and QCD sumrules W\ f B Jf B = 1.16 ± 0.05 and f Ds /f D = 
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Input parameters 







(220 ± 30) MeV 


J n 




93 MeV 


(aa) 




-fO 240 GeV) 3 


a (A ) 




0.50 






0.36 GeV 2 


f R fn,J2Mnci B A * 


/ J 7Y 


(0.64 ± 0.06) GeV 2 


ctairrih) 




0.21 


777 » 




4.8 GeV 






1.1 


C*(A V ) 

v V X/ 




0.05 






0.85 






fl 95 


u 2 (l>) 




0.30 GeV 2 


fn fn*V2Mna^' 




(0.51 ± 0.05) GeV 2 






0.36 


m c 




1.4 GeV 


q(A x ) 




0.9 


C, C (A X ) 




0.08 


C C M (A X ) 




1.15 


C C K (A X ) 




0.75 



TABLE I: Input parameters of HLxQM in the B- and D- sector. 

1.15 ± 0.04. From table [TT] and [TTT| , we see that we are a little low in the 5-sector but 
in the .D-sector the result agree nicely. The decay constant fn s have been measured 



f Ds = (264 ±15 ±33 ±2 ±4) MeV and P§ f Ds = (280 ± 19 ±28 ±34) MeV. It is somewhat 



higher than our result, but we are within la of the experimental result. The ratio fn*/ fn 
has been calculated in HQEFT sum rules with the result ||0|] — 1.07 ± 0.03 and 

Jd*/ Id — 1-37 ± 0.04, which also agrees perfect with our results. 

The coupling g^* Dn has been measured |4l| g®* Dn = 0-°9 i 0-01 ± 0.07. Our prediction 
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0.55 i— 

0.50 - 

0.45 - 

0.40 ^ 

0.35 - 

0.30 - 

0.25 - 
0.20 

0.15 - 
0.10 

0.05 



<)A 

9 a(B*Btt 
g A (B*B*Tr 
J L 



0.18 0.20 0.22 0.24 

m (GeV) 



0.26 



0.28 



0.30 



FIG. 8: The strong coupling constant in the 5-sector 



1.00 



0.80 - 



0.60 - 



0.40 - 



0.20 - 



0.00 



_L 



0.18 0.20 



9A 

9a(D*Dtt) 
g A (D*D*Tr) 
J l_ 



0.22 0.24 

m (GeV) 



0.26 



0.28 



0.30 



FIG. 9: The strong coupling constant in the -D-sector 



agrees well with this result. The experimental result has also been predicted by a bag model 
calculation g|, 9a ^ = °- 60 and 9%** = °- 57 - 



In conclusion, we have constructed a model which gives a reasonably good description 
of decay constants, the chiral axial coupling g A , masses and the Isgur- Wise function. We 
observe that the coupling g A (leading order and corrected) is smaller in the 5-sector than 
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> 

CD 

o 



•3 



0.10 
0.18 



i 1 " r~ 

Leading order (LO) 
LO +l/m Q 
LO + X + l/mq 




_L 



0.20 0.22 0.24 0.26 

m (GeV) 



0.28 



0.30 



FIG. 10: Jb as a function of m 




FIG. 11: as a function of m 



for the D-sector. Furthermore, we find that fjQ > \jf G both in the B- and the Z)-sector [43 



We have also showed that it is possible to systematically calculate the I/tuq corrections. 
For the decay constants in both the B- and D-sector they are of the same size as the chiral 
corrections as, can be seen from figure pT0|-p~3|. 

The model may be used to give predictions for other quantities. Especially, it will be 
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Predictions of HLxQM 



G B 


(7.7 ±0.6) GeV 1/1 




(0.315 ± 0.020) GeV 


91 


(1.3 ±0.2) GeV 


92 


(0.39 ± 0.05) GeV 


£\ 


-(0.7 ±0.2) GeV 


Ai 


1.0 ± 0.2 


2 


(0.41 ± 0.02) GeV^ 


9A 


(0.42 ± 0.06) 


B* Bit 

9 A 


(0.31 ± 0.11) 


B* B*TT 

9A 


(0.22 ± 0.13) 


Sb 


(165 ± 20) MeV 


Ib* 


(170±25) MeV 


fBs 


(170±20) MeV 


Jb* 


(175±25) MeV 


f B */fB 


l.06±0.03 


!bJ!b 


l.07±0.02 



TABLE II: Predictions of HLxQM in the B- sector, the errors in the predictions is a consequence 
of the error bars in the input parameters. 



suitable for calculation of the 5-parameter for B — B mixing [IT3 . 
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Predictions of HLx QM 



Gn 




(6.8 ±0.4) GeV -1 / 2 


/ a s 

\ 7T 




(n gnn _i_ n nnn) n P V 

\ \J •tJ\J\J —1— VJ . VJZj VJ ( VJC V 






(0.8 ±0.1) GeV 






(0 qq _i_ n rvA HpV 
i u.oo in u.u^ti ucv 


=•1 




i u.u in u.^ J vjc v 


Ai 




0.7 ± 0.1 


u 2 

r 7T 




(0 32 ± 03) GeV 2 






55 ± 08 

V / • J .J — 1 — V ' • V J v./ 




Dn 


0.46 ± 0.15 




D*TT 


(0.27 ± 0.22) 


fn 




C1Q0 ± 20) MeV 


p 

Id* 




(220 ± 35) MeV 


Id s 




(205 ±15) MeV 






(235 ±35) MeV 




IfD 


1.22 ±0.09 




IfD 


1.16 ±0.03 



TABLE III: Predictions of HLxQM in the D- sector, the errors in the predictions is a consequence 
of the error bars in the input parameters. 

APPENDIX A: LOOP INTEGRALS 

The divergent integrals entering in the bosonization of the HLxQM are defined : 

d d k 1 



h 



f d d k 1 

J (2ir) d k 2 -m 2 { } 

f d d k 1 

3/2 = J {2ix) d {v k){k 2 - m 2 ) ( } 

f d d k 1 

2 ~ J (2n) d (k 2 - m 2 ) 2 ( j 




FIG. 12: fs* as a function of m 




FIG. 13: fo* as a function of m 
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The defined as: 



m _ . m 2 



l07T 8ft 2 

_.8-3tt _.3tt-4 

Kl= *384m 3 K2=Z 384m 2 ^ ^ 

= , 8 + 3tt = ,2-3tt 

K3_Z 384m 3 K4_ *192m 2 ^ ^ 

"■^ss?- (A7) 

Integrals involving a heavy quark propagator and a light quark propagator can be sym- 
metried using the following formula : 

d d k 1 



in, A 



(2ir) d (k 2 - m 2 )p(v ■ k - A) c < 



r(p)r(g) 7 7 (2vr) d (k 2 - m 2 + 2Xv ■ k)P+i 
o 

There are three interesting limits, where this integral can be written in a rather compact 
form : A, m — > and A = m : 



,„,, , v ,_, , -i)^»r( g /2)r(p + g /2-d/2) /^y+«/ 2 - d / 2 
m ( 4vr )d/2 r(p)r(g) \™ 2 ' 



r o,A (-1)^ r(d/2- P )r(2 P + g-rf) / i y+^/ 2 
p>« - 2 2 P -d( 4vr )d/2 r(p)r( g ) V A2 " 



-lf +q i r(p + g/2 + 1/2 - d/2)r(p + g/2 - d/2) / 1 



p+q/2-d/2 



jm,m _ \ \i j. I I I J yr_ 1 ~ "7 | J_ (Alii 

M (47r) d / 2 r(p + g + l/2-d/2) \m 2 ) V J 

where d = 4 — 2e is the dimension of space. As a check equation ( |A9|) gets the well known 



is 



form in the limit q — > (T(q/2)/T(q) — > 2),L™q =J p . In the limit p — > 0, = 0, this i 
because there is no mass scale entering in the integral. 

In the general case for A,m ^ 0, there is no compact form of L™' q A , but all integrals 
needed in calculations can be relatet to the following integrals : 

L™[ A = ^(-- ln(m 2 ) + 2 - 2F(m/A)) (A12) 
87r V e I 
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Ag^ + BvW 



A 



(2ir) d (k 2 - m 2 )(v ■ k - A) 

1 f d d k k 2 -{vk) 2 



d-lj (2n) d (k 2 - m 2 )(v k - A) 



/A \-- + ln(m 2 ) - l)(m 2 - -A) - -F(m/A)(A 2 - m 2 ) - -(m 2 - -A 2 )Ml3) 



16tt 2 I v e v " ' /v " 3 7 3" v "' /v "" y 3 V " 6 

, d d k (vkf 
B = -A+ ' 



(27r) d (k 2 -m 2 )(vk - A) 

/(_! + l n ( m 2 } _ 1)(2m 2 _ 8 A) _ 4 F(m/A)(4A2 _ m2 , 
167T^ I e 3 3 

-l(m 2 - ^A 2 ) \ (A14) 



where: 



— y/x 2 — 1 tan 1 (\/x 2 — 1) x > 1 
= <( (A15) 

Vl - tanh"^^! - x 2 ) x < 1 



In the case of A > m we have ignored an analytic real part in ( |A12[ ). As a check on these 



calculations we can use equation ( |A11| ). Equation ( |A12| ) coincides with the one obtained in 
21| however equation (|A14 ) differs by a factor — 2/3(m 2 — 2/3A 2 ) inside the parenthesis of 



the expressions for A and B. This is presumably due to the factor l/(d — 1) = (1 — 2/3e)/3 
in A. 

In the case of the IW function we also need the following integral : 

, d d k 1 

L(p,u) 



(2n) d (v ■ k)(v' ■ k){k 2 -m 2 )P 
2i(-l) p T(p+l-d/2) 



, ,„r(uj) where 
(4n) d / 2 T{p)(m 2 Y+ l - d / 2 V ; 

r(uo) = - 1 In (u + Vu 2 - l) (A16) 
yuo 2 — 1 V ' 

APPENDIX B: SU(3) TRANSFORMATION PROPERTIES 

The Lagrangian of the light quark sector is 

£> X QM = Qdl ■ Dq L + q R i^ ■ Dq R - q L M q q R - q R M\ q L ~ m(q R T)q L + q L Eq R ) , (Bl) 
The left - and right - handed projections of the quark fields transform as 

Ql ~^V L q L , q R -> V R q R , (B2) 
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where V^e SU(3)l and V^e SU(3)r. The octet meson field £ transforms as 

S^V L SV^. (B3) 
The £ field transforms more complicated as 

e-i/£y£ =y^ (B4) 

where UeSU(3)y. The constituent quark fields xl — <^<?l an d Xr — £qri an d the heavy 
meson field transform in a simple way under SU(3)y: 

Xl^Uxl , Xr^Uxr , H V ^H V UK (B5) 

The vector and axial fields transform as 

V„ -> U V„ f/ f + zt/fy C/t , ^ _> [/^ f/t . ( B6 ) 

Note that the weak current in (p2[) transforms as 

^-^(Vj) (B7) 
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